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Abstract 

A methodology to improve the stability robustness of feedback control systems 
designed using direct eigenspace assignment techniques is presented. The method consists 
of considering the sensitivity of the minimum singular value of the return difference 
transfer matrix at the plant input to small changes in the desired closed— loop eigenvalues 
and the specified elements of the desired closed— loop eigenvectors. Closed— form expressions 
for the gradient of the minimum return difference singular value with respect to desired 
closed— loop eigenvalue and eigenvector parameters are derived. Closed— form expressions 
for the gradients of the control feedback gains with respect to the specified eigenspace 
parameters are obtained as an intermediate step. The use of the gradient information to 
improve the guaranteed gain and phase margins in eigenspace assignment based designs is 
demonstrated by application to an advanced fighter aircraft. 

Introduction 

A fundamental objective in the design of flight control systems is to change the 
transient response of the flight vehicle to a desirable one using feedback control. As 
discussed in Ref. [1], the so called direct eigenstructure (or eigenspace) assignment 
techniques are really well suited to designing feedback control systems to meet this 
objective. Various applications of the eigenspace assignment techniques have appeared in 

Controls Engineer, Member AIAA. 


1 



the literature in the recent past. Some of these are - design of pitch pointing flight control 
systems [2], design of flutter and gust load alleviation systems [3], and transient response 
shaping for flexible vehicles [4]. All these applications were multivariable in nature and 
demonstrated that direct eigenstructure assignment is a viable multi— input multi-output 
(MIMO) control system design technique. 

One of the drawbacks of direct eigenspace techniques, as compared to some other 
multivariable techniques — specially the Linear Quadratic Gaussian/Loop Transfer 
Recovery (LQG/LTR) approach [5], is that the synthesis procedure does not guarantee 
stability robustness with respect to variations in plant dynamics. Even using the Linear 
Quadratic Regulator (LQR) based methodology to asymptotically approach the desired 
eigenspace [6] does not guarantee the well known stability margins of LQ Regulators [7] as 
the procedure results in non— diagonal control weighting matrices which violates the 
conditions under which the LQ Regulator stability margins are guaranteed. Also, as 
discussed in Ref. [4], the direct eigenspace assignment techniques are preferable to the LQR 
based approach because the LQR approach requires very high actuator bandwidths for the 
desired eigenstructure to be achieved. 

In direct eigenspace assignment techniques the design parameters are the desired 
closed— loop eigenvalues and specified elements of the closed— loop eigenvectors. Once the 
design parameters are specified, the feedback control gains are uniquely determined 
(provided enough parameters are specified — see Ref. [8] for discussion of limits on 
achievable eigenspace using direct eigenspace assignment). So, given a set of specifications, 
the feedback control gains will provide the desired closed— loop transient response (or come 
as close to it as possible within the system constraints), but they might result in a system 
with poor stability robustness, i.e. a small change in the plant dynamics may cause the 
closed— loop system to go unstable. The designer is then faced with the dilemma of how to 
change the design specifications such that the resulting feedback system will also provide 
adequate stability robustness. Note that in general the designer does have a certain amount 
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of freedom in choosing the design specifications — rarely does he want an exact value for a 
closed— loop eigenvalue or exact shape for a corresponding eigenvector, the specifications 
are rather in terms of desired regions for the closed-loop eigenvalues and acceptable sets of 
eigenvector shapes. The objective of this paper then is to develop a methodology which will 
provide adequate information to the designer to change the design specifications in a 
systematic step-wise manner such that at each step the guaranteed stability robustness of 
the feedback system is improved while the eigenstructure is within the desirable regions. 

In multivariable feedback systems, a reliable (but sometimes conservative [7]) 
measure of stability robustness is the minimum singular value of the return difference 
matrix at the plant input evaluated as a function of frequency. The methodology presented 
in this paper is based on sensitivities of the minimum singular value of the return difference 
matrix to the design parameters, which in this case are the desired closed— loop eigenvalues 
and eigenvectors. Note that the notion of using return difference singular value sensitivities 
to design robust controllers is not new. Singular value sensitivities to compensator 
parameters were used in Ref. [9] to directly design robust reduced— order compensators, and 
singular value sensitivities to plant parameters (elements of the plant system matrices) 
were used in Ref. [10] to determine which elements need to be modeled more accurately for 
the feedback system to guarantee stability. 

In the following a technique for solving the feedback gains for direct eigenspace 
assignment is first briefly described and the mathematical problem formulation for deriving 
analytical expressions for the return difference singular value sensitivities is presented. 
Closed— form expressions for the singular value gradients with respect to closed— loop 
eigenvalues and elements of the closed-loop eigenvectors are then developed by first 
considering the singular value gradients w.r.t. the feedback gains and then deriving the 
closed— form expressions for the gradients of the feedback gains to the closed— loop 
eigenspace parameters. Finally, the use of the gradient information to improve guaranteed 
stability robustness is demonstrated by application to a modern fighter aircraft. 
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After the completion of the present manuscript, it was brought to the author's 
notice that somewhat similar research has been reported earlier in Ref. [11]. However, as 
pointed out in a later section, the research reported herein differs from that reported in 
Ref. [11] in the important detail of the choice of independent parameters with respect to 
which the singular value sensitivities are calculated. 

Problem Formulation 

Direct Eigenspace Assignment Gain Synthesis 

In the direct eigenspace assignment technique, the control objectives are stated in 
terms of a desired eigenstructure for the augmented system. For the full— state feedback 
case, the synthesis problem is as follows : 

Given a linear, time— invariant dynamical system with the state— space 
representation 

x = Ax + Bu (1) 

with xdR n and udR m , find a control law of the form 

u = -Kx (2) 

to achieve some desired eigenspace for the augmented system 

x = (A-BK)x (3) 

To determine the feedback gains K, note that the augmented (closed— loop) system 
eigenvalues and eigenvectors are related by 

(A— BK)v = AjV , i = l,...n (4) 

where A. is the i^ closed— loop eigenvalue and v is the corresponding closed— loop 
i c i 

eigenvector. 

For full— state feedback, the limitation on the achievable eigenspace is that all the 
desired closed— loop eigenvalues can be exactly placed while only "m" elements of their 
associated eigenvectors can be exactly achieved [8] (here, m = dimension of u and it is 
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assumed that the system given by (1) is controllable). Since in general m<n, we cannot 
exactly obtain all elements of the desired eigenvector for each closed-loop mode. One 
approach for determining the feedback gains is to obtain the "best" achievable 
eigenvectors, for each of the closed— loop modes, so as to minimize the mode's cost function 
Jj given by 



i _ _ * 

' K -V H ) 


Qi < v a- v d.) 
1 1 


i = l,...n 


( 5 ) 


where 


— fin 

v = r n achievable eigenvector associated with eigenvalue A- 

dj ♦ 1 

til 

v d = i in desired eigenvector 
th. 

Qj = i 11 n— by— n symmetric positive semi-definite weighting matrix on 
eigenvector error elements, 

* 

and [•] denotes complex— conjugate transpose of [•]. 

Equation (4) can be rewritten as 


(A.I-A)v = " BKv a. ( 6 ) 

Defining the vector w. = — Kv and using eqn. (6), the solution to minimizing the cost in 

i a j 

eqn. (5) is obtained as (see Ref. [4] for a complete derivation) : 

= [L*QjLj] _ 1 L*Q i v di (7) 

where L = (AjI-A) _1 B. 

Once Wj are obtained, the achievable eigenvectors are given by 

v a = LjWj , i = l,...n (8) 

and the feedback gains are obtained as 


K = -WV -1 

where W = [w, w 9 .... w_l and V = [v„ v 0 .... v 0 ]. 

1 Z D. 


( 9 ) 


Note that this algorithm requires that the specified closed-loop eigenvalues Aj be 
distinct and different from the open-loop eigenvalues (eigenvalues of plant system matrix 
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A). Also, we are considering full-state feedback rather than reduced-order output feedback 
because reduced— order output feedback further restricts the achievable eigenspace [S], and 
full-state feedback can always be implemented with state estimation without any 
significant loss in stability robustness by using either the loop recovery procedure of the 
LQG/LTR approach or an eigenspace assignment based loop recovery procedure discussed 
in Ref. [3] . 

Problem Statement 

Given a multivariable system as in eqn. (1) with a control law of the form (2) and 
the state feedback gains given by eqn. (9) such that the closed— loop eigenvalues are 


A 2i-1 2i = “V-^i 1 = lv " k ’ ^ and ^ >0 

and Aj = — ^ , i = 2k+l,...n (10) 

and the desired eigenvectors are 

\i-l 2i = fof A— l,2i = -V^l 

and v d = ool.fi/j-] for \ = ( n ) 

u i J 

we wish to derive analytical expressions for the sensitivities of the minimum singular value 

of the return difference matrix at the plant input to the specified eigenstructure elements, 

i.e. closed— form expressions for the partial derivatives 

<9u[I+KG(s)] 


( 12 ) 


where G(s)=(sl— A) 1 B with s being the Laplace operator, a[-] denotes minimum singular 
value of [•], and £ represents the eigenspace parameters 7., £., 77. , p ^ and i/- defined in 


(10) and (11). 

In Ref. [11], the parameters with respect to which the singular value sensitivities are 
calculated are the desired closed— loop eigenvalues Aj (and hence 7^ and ^ as above), and 
the elements of the vectors w^ given by Eqn. 7. (wj corresponds to the t^ of Ref. [11]). As 
seen from Eqn. 7, wl depend on the choice of A- , and hence it is inappropriate to choose w^ 
and Aj as independent parameters for calculating the return difference singular value 
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sensitivities. It is important to note that the design specifications available to the control 

designer are in terms of the desired eigenvalues and eigenvectors, A- and v , respectively, 

a. 

and not in terms of Aj and Wj (tj). Therefore, it is more meaningful to choose Aj and v^ as 

the independent parameters with respect to which the stability robustness singular value 
sensitivities are obtained. 

The discussion in the present paper is limited to singular value sensitivity analysis 

for the return difference matrix at the plant input. In general, it is important to guarantee 

the stability robustness at the plant output also. As will be apparent from later discussions, 

the ideas presented in this paper can be extended to deriving sensitivities of the minimum 

singular value of the return difference matrix at the output simply by replacing I+KG(s) 

with I+G(s)K. Also note that the singular value sensitivities to the natural frequency and 

damping, u^ and ( respectively, of a complex closed— loop mode can be obtained by using 

da_ da <97j da dd j 

du dj- du d6- du 

n i 1 n i 1 n i 

and 


da da <9 7. da dd^ 

dCj H x dS. dq 


with 




du 


n i 


dr 

= Cj ; — 

a Cj 


= u 


n i 


d 6 . , dd, ^i^n. 

I a a 


du 


n- 


8c i Jwf 


where we have used the equalities 7. = C-u , and 6 - = u J 1— C?. 

'1 u 1 nj 1 b i 

Lehtomaki et al. [ 7 ] have shown that if 

a[I+I<G(;h;)] > a Q , 0 <u <00 ( 13 ) 

for some constant a Q < 1, then simultaneously in each loop of the feedback system there is 
guaranteed gain margin (GM) given by 
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( 14 ) 


GM = — 
l± a o 

and also a guaranteed phase margin (PM) given by 


PM = ±cos 



2 


(15) 


Therefore the gradient information provided by expressions of the form (12) can be used to 


move the eigenspace parameters within the desirable regions in a systematic, step-wise 
manner such that the lowest value of the minimum singular value of the return difference 
matrix, and hence the guaranteed gain and phase margins, are improved at each step. 


Singular Value Sensitivity Derivation 

In Ref. [12] it has been shown that for a general complex matrix H of rank m which 
has distinct singular values <r-, i = l,...m, the sensitivity of the singular value <j. with 
respect to a real parameter p is given by 

^j[H] 


dp 


= Real 


-* <9H - 

u i — v i 
. ' dp 


( 16 ) 


where w and u^ are the right and left singular vectors [13], respectively, corresponding to 
the singular value 

The eigenspace parameters defined in (10) and (11) are all real and the open— loop 
state frequency response matrix (G(jui)) is independent of these closed— loop eigenspace 


= Real 


(17) 


parameters. Therefore, using (16), we get 

<9a[I+KG(jw)] 

where v and u are the right and left singular vectors, respectively, of [I+KG(jw)] 
corresponding to the minimum singular value a. With I< given by (9), we have 


u — G(ju) v 
. H 


dK 


m t K dv 


V 


-1 


(18) 


[dt d(\ 

The expressions for the feedback gain sensitivities to eigenspace parameters are fully 
expanded in the Appendix by first considering the eigenvalue parameters (r/., 7^ and ^) and 
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then the eigenvector parameters (iaj, ^ and /9y). An algorithm which uses the sensitivity 
information for stability robustness improvement is briefly discussed in the following 
section and finally an example application is presented. 

Stability Robustness Improvement Algorithm 

Based on the derivation of the return difference singular value sensitivities to 
eigenspace parameters in the previous section and the Appendix, a step— by— step algorithm 
for improving the guaranteed stability robustness of direct eigenspace assignment feedback 
control designs is as follows : 

Step 1 : Formulate the eigenspace requirements; the desired closed— loop eigenvalues 

and closed-loop eigenvectors, and identify the design freedom available in 
the closed-loop eigenspace specification. 

Step 2 : Solve for the control feedback gains K, and check for the guaranteed stability 

robustness using the minimum singular value of the return difference matrix 
at the input, a[I+KG(jw)]. If the guaranteed stability margins are 
acceptable, then stop, otherwise go to step 3. 

Ot t r At r 

Step 3 : Calculate — — and — for the eigenspace parameters £, identified in Step 1, 

n ds 

for which some design freedom is available. The calculation procedure is as 

discussed in the Appendix for each case of eigenspace parameter. 

8K 0£[I+KG(jw)] 

Step 4 : Using the results of Step 3, calculate — from (18) and then 

di 

from (17). 

Step 5 : Using the information from Step 4, make small changes in the eigenspace 

parameters such that a[I+KG(jw)] will increase in the desired frequency 
region while making sure that the changed eigenspace parameters are within 
the bounds specified in Step 1. Go to Step 2. 

The above algorithm can be implemented in the form of a constrained optimization 
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technique with the objective of minimizing the area under a specified value for the 
minimum return difference singular value frequency response within the constraint that the 
design parameters stay within certain specified bounds. Such an optimization approach is 
discussed in Ref. [9]. However, the choice of changing the eigenspace parameters in step 5 
requires a considerable amount of "engineering intuition" which cannot be easily put in the 
form of a computer program. Therefore, for the present study, a computer program was 
developed just to implement steps 2 to 4 and the selection in step 5 was based on 
knowledge about the particular dynamics being controlled. 

Example 

The flight vehicle model considered is the short period approximation for the 
AFTI/F— 16 aircraft as discussed in Ref. [2]. The model is for a flight condition 
corresponding to an altitude h = 3000 ft and Mach number M = 0.6. The equations of 
motion are given in the form of (1) with 

x = [t, q, «, <$ e , 

where 

7 = flight path angle 
q = pitch rate 
a = angle of attack 
6 q = elevator deflection 
<!>£ = flaperon deflection, 
and 

“ = [*«• y T 

where 

6 = elevator deflection command 

CL 

S { = flaperon deflection command. 

The system matrix A and the control distribution matrix B have the following numerical 
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values : 


'0 0.00665 1.3411 0.16897 0.25183 ] [0 0' 

0 -0.86939 43.223 -17.251 -1.5766 0 0 

A = 0 0.99335 -1.3411 -0.16897 -0.25183 ; B = 0 0 

0 0 0 -20 0 20 0 

0 0 0 0 -20 J 0 20 

The eigenvalues of the open— loop system are given by 
A t =-7.662 ] 

• unstable short period mode 

A 9 = 5.452 J 

Ag = 0.0 pitch attitude mode 
A^ = —20 elevator actuator mode 
Ag = —20 flaperon actuator mode 

The control design objective is to provide decoupled tracking of flight path and 
pitch attitude commands with a well damped response and zero steady— state error to step 
commands. The control law is of the form 


u = Fy c -Kx (19) 

with y c = [ 7 c , # c ] . F is the feed— forward gain matrix and K is the feedback gain matrix. 
The feedback gains are to be obtained using direct eigenspace assignment techniques with 
the objective of providing decoupled flight path and pitch response modes. The feedback 
gains should also be such as to guarantee gain margins of at least ±3.5 dB for simultaneous 
gain changes in each control loop at the plant input and phase margins of at least ±30 deg 
for simultaneous phase changes in each control loop. These stability margin specifications 
translate into the requirement that ajl+KGQo;)] > 0.51. Once a set of feedback gains that 
satisfy the transient response and stability robustness requirements are obtained, the 
feedforward gains (F) will then be obtained using a special case of Broussard's command 
generator tracker [14]. As discussed in Ref. [2], a set of feedforward gains that provide 
command following with zero steady— state error is given by 
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1 


( 21 ) 


n = 

Ai 

fi 12 



A B' 


fi 21 

^22 


H 0. 


where H is defined by 

y = Hx (22) 

with y being the controlled variables of the plant. For the present example, y = [7, v\ ■ 

A set of feedback gains was obtained initially to achieve the desired eigenspace 
listed in Table 1. In Table 1, the short period frequency and damping (corresponding to 
A 1 9 ) were chosen to be w -7 rad/s and ( = 0.8, respectively, so as to meet the 
MIL— F— 8785C [15] specifications for Category A, Level I flight, and the flight path mode 
(Ag) was chosen to provide adequate bandwidth for flight path control. The eigenvectors 
corresponding to the short period mode and the flight path mode were chosen to minimize 
the coupling between the pitch rate and the flight path angle. The actuator mode 
eigenvalues were chosen to be close to their open-loop values and the corresponding 
eigenvectors were chosen to minimize actuator cross— feed. Note that we are only specifying 
two elements for each desired eigenvector ('X' in Table 1 denotes arbitrary), and since we 
have two control inputs the specified eigenspace can be exactly achieved as seen from the 
achieved eigenvectors listed in Table 1. The control feedback gains for this initial design 
are also listed in Table 1. The minimum singular value of the return difference frequency 
response at the plant input is shown in Fig. 1. We note from Fig. 1 that the minimum 
singular value at low frequency is much lower than that required to guarantee the desired 
stability margins. Therefore, although this initial design will meet the performance 
requirements, a control law redesign is necessary for the stability robustness specifications 
to be met. 

The singular value sensitivity procedure was then applied to the initial feedback 
design. The eigenspace parameters for which there is design freedom are : 

(a) Short period frequency and damping (and hence A^ 9 ). Based on MIL— F— 8785C, the 

allowable regions for Level I flight response, at the given flight condition, are 0.35 < 
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£ < 1.35 and 2.5 < o» n < 8.5 rads/sec. 

(b) Flight path mode (A^). Any value > 1 rad/s will provide adequate flight path 
control, however there will be an upper limit due to maximum allowable flap 
deflection. 

(c) The desired eigenvector elements v, (and corresponding element v, ) and v , , 

U 11 a 21 a 32 

i.e. the short period contribution to the flight path and the flight path mode 
contribution to pitch rate. Although these elements were chosen to be 0 (zero) in 
the initial design, the only requirement is that they be "small" (<< 1) in order to 
keep the coupling of the modes to a low level. 

There is no design freedom available in placing the actuator modes as it is desirable to keep 
these close to their open— loop values. 

The predicted change in the minimum singular value of the return difference matrix 

for a 10 % (percent) increase in the short period frequency and damping, obtained using 
da da 

— and — respectively, is shown in Fig. 2 and that for a 10 % increase in the flight path 
du n d( 

mode (| A^|) is shown in Fig. 3. The return difference singular value sensitivities to the real 
and imaginary parts of the desired eigenvector element v^ , and p ^ respectively, are 

shown in Fig. 4 and the sensitivity to eigenvector element v^ (j^ 9 ) is shown in Fig. 5. 

32 

From these figures we note that changing the eigenvector element v , , decreasing m 1 and 

a n 11 

increasing will be most effective in increasing the minimum singular value of the 
return difference matrix. 

Based on the results in Fig. 5, letting v, = -0.01+/). 01, and keeping the rest of 

a ll 

the desired eigenspace parameters the same as in Table 1, should result in an increase of 
sO.115 in the lowest value of the minimum singular value of the return difference matrix 
while still keeping the coupling between flight path and pitch rate, for the short period 
mode, to be very low. The minimum return difference singular value with the feedback 
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gains corresponding to this choice of eigenspace parameters, referred to as trial 1 from here 

onwards, is shown in Fig.6. Also shown in Fig. 6 is the minimum singular value predicted 

da da 

using the singular value sensitivities, i.e. a[I+K G(^y)] — 0.01 1- 0.01 — where K Q 

dfi 1{ dp n 

are the feedback gains corresponding to the initial design. Note that although there is not 

very good agreement between the actual and predicted values at low frequencies, the 

sensitivities did accurately predict the direction of the change (increase) in the singular 

value. Also note that the actual feedback gains for trial 1 and those predicted using 

— and — , both listed in Table 2, are virtually identical. From Fig. 6 we get 
dp n dp n 

a[I+KG(jw)] > 0.4 for trial 1, which is much improved over the initial design but still not 
high enough to meet the stated stability robustness requirement. Therefore we need to 
further change the eigenspace design parameters in order to improve the guaranteed 
stability margins. 

For trial 1 the singular value sensitivity calculations showed that the most effective 

way to increase the minimum singular value of the return difference matrix was to further 

decrease and increase p However, doing so will result in increased coupling in the 

flight path and pitch rate response for the short period mode which will be undesirable. 

Next to the desired eigenvector element v- , the return difference singular value was most 

Q 11 

sensitive to changes in the short period mode. The changes in a[I+KG(ju>)] from that for 
trial 1, using singular value sensitivities for a 10 % increase in short period frequency and a 
10 % increase in damping are shown in Fig. 7. From Fig. 7 we note that increasing both the 
short period frequency and damping by 10 % each over that for trial 1 will result in an 
increase of 0.114 in the minimum singular value of the return difference matrix. We will 
then have <r[I+KG(ju;)] > 0.51 which will satisfy the design requirement. 

The short period frequency and damping for a 10 % increase over that for trial 1 are 
w n = 7.7 rad/s and ( = 0.88 which are within the region for Level I flight requirements. 
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The desired and the achieved eigenspace for this change are summarized in Table 3 and the 
corresponding control feedback gains are also listed there. The resulting minimum singular 
value of the return difference matrix is shown in Fig. 8. From Fig. 8 we note good 
agreement between the actual singular value response and that predicted using singular 
value sensitivities. For this case also there was excellent agreement between the actual 
feedback gains and those predicted using feedback gain sensitivities for trial 1. The plot in 
Fig. 8 shows that this set of feedback gains satisfy the stability robustness design 
requirement, therefore this case will be referred to as the final design. The closed— loop 
performance of the initial and the final designs is compared in the following. 

The feed— forward gains for the initial and final designs, calculated using eqns. (37) 
and (38), are listed in Table 4. The response of the initial design to a unit step gamma 
command (7 (t) = 1 deg) is shown in Fig. 9. Shown in Fig. 9 are the time histories of the 
controlled variables 7 and 9 , and also the control input deflections £ and 6 ^. The response 
of the final design to a unit step gamma command was identical to that of the initial 
design, so those time histories are not shown here. From Fig. 9, then, we note that both the 
designs provide well-damped tracking of flight path commands with a reasonably fast rise 
time and without any perturbations in the pitch attitude. The responses of the initial and 
final designs to a unit step pitch attitude command (9 (t) = 1 deg) are shown in Figs. 10 

V 

and 11 respectively. Note that both the designs provide fast tracking of pitch attitude 
commands, however the initial design does so with no perturbation in the flight path angle 
while the final design does result in a small initial perturbation in the flight path angle.. 
This coupling of the flight path with the pitch attitude is due to small contribution of the 
short period mode to the flight path that was allowed in order to improve the guaranteed 
stability robustness. Therefore, in just two iterations the judicious use of the singular value 
gradient information led to a feedback control law design with much improved stability 
robustness while still maintaining acceptable performance. 
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Summary and Conclusion 

A methodology to improve the stability robustness of feedback control systems 
designed using direct eigenspace assignment techniques was presented. A full— state 
feedback gain synthesis technique for direct eigenspace assignment was briefly discussed 
and closed-form expressions for the sensitivity of the minimum singular value of the return 
difference matrix, at plant input, to changes in desired closed-loop eigenvalues and 
specified elements of the desired closed— loop eigenvectors were derived. The closed— form 
expressions for the sensitivity of the control feedback gains were obtained as an 
intermediate step. An algorithm discussing the steps involved in calculating the return 
difference singular value gradients was presented and the use of the gradient information to 
improve the guaranteed gain and phase margins was demonstrated by application to an 
advanced fighter aircraft. The aircraft example consisted of the short period approximation 
of the longitudinal dynamics with the design objective of providing decoupled tracking of 
flight path and pitch attitude commands using the elevator and flaperon as control 
effectors. The nominal feedback control design which provides ideal decoupling between 
flight path and pitch rate was found to have very low guaranteed stability robustness. 
Using the singular value gradient information it was shown that by allowing the short 
period mode contribution to the flight path angle to be non— zero, but still small (<< 
contribution to pitch rate), and by increasing the short period mode frequency and 
damping by 10 percent from that for the nominal design, the guaranteed stability margins 
could be increased significantly while still maintaining acceptable performance. 

Appendix 

Feedback Gain Sensitivity to Eigenvalue Parameters 
SK 

— : We first consider the n— 2 k real closed— loop eigenvalues, — 77. , i = 2 k+l,...n. Then 

‘ 
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dW 

drj. 


o| |o 

i d V l i . 

i^col. 


(A.l) 


because, as seen from (7), Wj, j^i, does not depend on Vy Using (7) we further have 

dWj _ ^ L T Q i L j J 1 t T^~ , rrT, 

dv j dr] { 


L i Q i 7 d. + t L i Q i L i 


1 5lT 

- 1 1 Q ; F 


n id- 

dr] j 1 


(A.2) 


which can be expanded to give 


3n , ' lW "' 


a? _ a. 

— - Q-L. + L • Q. — 
or] j dv- x . 


[L r Q.L.] 1 L T + — 

L 1 ^1 1 J 1 o 

dr] j 


a T 


Q;V 


i v d. 


Furthermore, using the definition of L-, we have 

<9L. 


1 = (-Vjl-A) 1 L i 


5?/. 


(A.3) 


(A. 4) 


d\v, 


Substituting (A. 4) into (A.3) we will get an expression for — i in terms of known 


dr]. 


quantities, and substituting the result in (A.l) we will get 


Next we determine 


dv 

dr] j 


3W 

dv j 


where V is the matrix of achievable eigenvectors as defined 


earlier. Noting that v , j^i, (as defined in (8)), does not depend on 7 /., we have 

j 1 


dV 

dv j 


' 1 dF 1 

0 1 — li l 0 


L I 


1 ! 


(A.5) 


and 


dv dL ■ «— 

— -i = — - w- + L. — i 

dv j dv^ dv- x 


dL, 


(A. 6) 




_ — , Cl v » 

Knowing — from (A. 4) and — i from (A.3), we can determine — -i from (A. 6) and then 

dV i dv [ dv { 


Or r 

substitute in (A.5) to get — . Once and — are known, — are determined from eqn. 


m 


ev 


dK 


dv, 


dv j dv j 


dv i 
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( 18 ). 

— : Next we consider the k complex modes, — 7j ± $p i=l,...k, by first deriving the 

feedback gain sensitivities to the real part of the complex eigenvalues. If we arrange the 


columns of W such that 


W = [wj, w 2 , ... w 2i _ p w 2i , ... w 2M , w 21 , Wj (j=2k+l,...n)] 
where w 2 \—\ corresponds to ^2i-l = ~ an d w 2i corresponds to A^= the 11 

w 2i = conj(w 2i _ 1 ) where conj(-) denotes complex conjugate of (•) and we have 


dw 

9 7 j 


0 ! — 2i— 1 ! — 2i I 0 

■57. 1 dj- 1 

1 1 1 1 1 


(A.7) 


ottt 

Proceeding just as in the case of above, we get 

d r] j 


^2i 1 - [ L 2i-iQ2i-l L 2i-l^ 


^7j 


^2i-l 


dj- 


*T* 

^2i-l L 2i-l +L 2i-1^2i-l 


gL 2i-r 

^7j . 


* 

* 1 * ^2 i -1 1 

t L 2i— 1 ^ 2i— 1 L 2i— 1 ^ L 2i— 1 + 


and 


5L o: 


2 i -1 


07i J 


= {(-7^)1 - A} 1 L 2i _, 


Q 2 i -Ud 2i _ 1 


(A.S) 


d 7i 


(A.9) 


g- 5[conj(w 2i _ 1 )] 

Using the relationship — 2i = = conj 

j 

use of (A.9), (A.S) and (A.7). 

Next, just as in (A.7), we have 


— 2i— 1 

l a 7j J 


dW 


, we can get - — by making 

d 7i 


OV 

dl [ 


1 dv 1 dv 1 
0 1 — =21—1 1 — 2i 1 0 


L l 


1 d 7; 1 dj . 1 

1 1 1 1 J 


(A.10) 


Also v = conj(v ) and from the definition of v in (8) we have 


a 2i 


2i — 1 
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S-i— j + L 9 j_j — 2i— 1 
»\ 2 ‘ 1 2 ‘ 1 »7j 

dW 


(A. 11) 


Using (A. 9), (A. 8) and (A. 11), we can determine — and hence determine an expression 


d 7; 


for — 


^7j 


— : The procedure for determining the control feedback gain sensitivities to the 

as i 


imaginary part of the complex closed-loop eigenvalues is similar to that for determinin 

— with y- replaced by 6 -. The only difference is in calculating — , which is now 

d'y^ 1 1 86, 


— ^ - A)} l = -r - 11 

86 . 1 1 


(A. 12) 


Using (A. 12) , instead of (A.9), in (A.8) and (A.ll), we can get analytical 

dv\ 


expressions for — 2i-l and — — 2i— 1 , respectively, and then we can determine — using 


8 K 


86 . 86 . 


86 .. 


(A. 7) and (A. 10) with 7j replaced by &. 

Feedback Gain Sensitivity to Eigenvector Parameters 

: We first consider the feedback gain sensitivities to the eigenvector elements v-- 

a "y 

. .... 8W , 8V 

corresponding to the real eigenvalues —77. for i = 2k+l,...n. The sensitivities and — 

1 d v ■■ dv-- 

ij 

are as given by eqns. (A.l) and (A. 5), respectively, with replaced by Vy. Furthermore, 
noting that U does not depend on elements of the desired closed— loop eigenvectors, we 
have 

(O' 


= [Uq^tUq. 

Sr.. 


it 

■ j n element 


(A.13) 


and 
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(A. 14) 


dv 


a -i = L. 


5v ij ‘ dv U 

The expressions (A. 13) and (A. 14) can be used to obtain - — and — , and then — can 

dv.. dv-. dv.. 

ij ij ij 

be obtained using (18). 

8K 

— : Next we consider the eigenvector elements corresponding to the k complex modes, 

% 

—'y.ijS., i=l,..k, by first deriving the feedback gain sensitivities to the real part of the 

“ij) 

respectively, with replaced by p,y j. Furthermore, we have 


eigenvectors (/!••). The forms for and — are obtained from (A. 7) and (A. 10), 

■'"ij a "u 


fa- 1 “ IbQihr'b^i 


a. 

■ j n element 


(A.15) 


and 


— ^2i— 1 = L„. . — 2i— 1 
3 "ij 8 "ij 


ry 

Using (A.15) and (A. 16) along with the equality — 2i = conj 

dp 


•J 


— 2i— 1 

l^n J 


(A.16) 
(and similarly for 


dv Hr, r A, r /Vf-a 

— — 2i), we can obtain — — and — , and then determine — . 

^ij d d V] dii.. dii.. 


dK 

9P ij 


ij 


: The procedure for determining the feedback gain sensitivities to the imaginary part 


of the desired complex closed— loop eigenvectors is similar to that for determining the 
sensitivities to the real part with //jj replaced by p-^. The only difference is in calculating 


d w 
dp 


■2i— 1, which is now 

ij 


f 2i -i = Mihr 1 ^ 

df> ij 


■ j^ 1 element 


(A.17) 
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Using (A. 17), instead of (A. 15), and proceeding as for the previous case, we can get 


and^ 


Or/ 

and hence obtain — 



dW 
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Table 1 Initial Design Eigensoace Summary 


State : x = [7, q, a, S Q , 6^ 


Eigenvalues : A^ 9 = — 5.6±j4.2, Ag = —1.0, A^ = -19.0, 


A. = -19.5 
0 


Eigenvectors : 
Desired 


Achieved 


% 

1 

il 

7 d3 ' 

'd. 

V , 

t dg 

V 

&1 

V3 

V 

S-4 

v„ 

as 


'0' 


V 


X 


X 



0 


1 


"— 0 . 006" 


'-0.014' 


1 


0 


X 


X 



1 


0 


1.072 


0.060 


X 


X 


X 


X 



-0.114tjQ.086 


-1 


-0.051 


0.011 


X 


X 


1 


0 



— 0.070Tj0.533 


-2.801 


1 


0 


X 


X 


0 


1 



0.629±jQ.S14 


3.234 


0 

L -1 


1 


Feedback Gains : K = 


-3.250 -0.891 -7.112 0.526 0.0S4' 
6.101 0.898 10.02 -0.420 -0.102 


Table 2 Feedback Gains for Trial 1 


Actual : I< = 


-2.865-0.827 -6.624 0.485 0.081 
1.833 0.196 4.612 0.035 -0.061 


Predicted : 


K = 


-2.865 -0.827 -6.624 0.485 0.081 
1.S30 0.195 4.609 0.035 -0.061 


23 


Table 3 Final Design Eigensoace Summary 


State : x = [7, q, a, 6 , <5f ' 


Eigenvalues: A^ 9 = — 6.78±; 3.66, Ag = — 1.0, A^ = —19.0, Ag = — 19.5 


Eigenvectors : 

Desired Achieved 


Vi 

d 1 ,2 

nEil 
1> 1 

V , 

d. 

Id = 

i 


v a 3 

cl»4 

% 

'-o.oi+jo.or 


rr 


X' 


X' 



■-0.010±j0.010' 


1 


'-0.006' 


— 0 . 014 

1 


0 


X 


X 



1 


0 


1.072 


0.060 

X 


X 


X 


X 



-0.104tj0.072 


-1 


-0.051 


0.011 

X 


X 


1 


0 



0.023tj0.414 


-2.801 


1 


0 

X 


X 


0 

1 

1 



0.637±j0.245 


3.234 


0 


1 


Feedback Gains : K 


'-3.280 -0.954-7.289 0.585 0.090 
0.138 -0.101 2.423 0.232 -0.043 


Table 4 Feedforward Gains 


Initial Design : 
Final Design : 


F = 
F = 


-0.373 -2.877' 
4.124 1.976 

— 0.373 -2.906' 
4.122 -3.984 
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MAGNITUDE 



FIGURE 3. - CHANGE IN £[ I+KGC j (J) ). FOR 10 PERCENT IN- 
CREASE IN |X 3 I (FLIGHT PATH MODE). 



FIGURE 4. - MINIMUM RETURN DIFFERENCE SINGULAR VALUE 
SENSITIVITY TO REAL AND IMAGINARY PARTS OF THE DE- 
SIRED FLIGHT PATH ELEMENT OF THE SHORT PERIOD MODE 


/do do \ 

EIGENVECTOR r— AND • 


MAGNITUDE 



FIGURE 5. - MINIMUM RETURN DIFFERENCE SINGULAR VALUE 
SENSITIVITY TO DESIRED PITCH RATE ELEMENT OF THE 

/ do \ 

FLIGHT PATH MODE EIGNEVECTOR I — I. 



FIGURE 6. - JJ[ I+KG( ju) I FOR TRAIL 1 - ACTUAL AND PRE- 
DICTED (USING GRADIENT INFORMATION). 
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0 .5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 

TIME, SEC 

FIGURE 9. - INITIAL DESIGN RESPONSE TO STEP y r = 1 DEG. 



0 .5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 



TIME, sec 

FIGURE 10. - INITIAL DESIGN RESPONSE TO STEP 0 C = 1 DEG. 
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6. (DEG) e, y (DEG) 



0 .5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 

TIME, sec 

FIGURE 11. - FINAL DESIGN RESPONSE TO STEP 9 C = 1 DEG. 
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